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Abstract 

We report a simple analytic solution for the velocity u of the transverse flow of QGP at 
a hadronization front in Bjorken's scenario. We establish scaling properties of the trans- 
verse flow as a function of the expansion time. We present simple scaling formula for the 
expansion velocity distribution. 
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Landau's hydrodynamic stage is a part of all scenarios for the evolution of the hot 
and dense matter (quark-gluon plasma - QGP) formed in ultrarelativistic heavy ion colli- 
sions . It is now well understood that Landau's complete stopping of Lorentz-contracted 
colliding nuclei is not feasible because of the Landau-Pomeranchuck-Migdal (LPM) effect, 
i.e., the finite proper formation time Tq (0, for the modern modern QCD approach to 
the LPM effect see [^, the early works on LPM phenomenology of nuclear collisions are 
reviewed in p), although evaluations of tq and of the initial energy density emax remain 
controversial [@]. The corollary of the LPM effect in conjunction with the approximate 
central rapidity plateau is the rapidity-boost invariance of initial conditions. The cor- 
responding solution for a longitudinal expansion in an 1+1-dimensional approximation, 
neglecting the transverse flow, was found by Bjorken ([0, see also [0). There is some 
experimental evidence H, |^, although a disputed one |n|], for a transverse flow which 
must develop if the lifetime of the hydrodynamical stage is sufficiently long. 

In this communication we present a simple solution of the Euler-Landau equation for the 
velocity of transverse expansion, u, gained in the hydrodynamic expansion of QGP before 
the hadronization phase transition. Our solution shows that for the usually considered 
lifetime tb of QGP the transverse flow is non-relativistic. It is only marginally sensitive to 
properties of the hot stage and offers a reliable determination of tb if the radial profile of 
the initial energy density is known. We find that the w-distribution is a scaling function of 
u/um, where maximal velocity of expansion. 

We start with the familiar Landau relativistic hydrodynamics equations 

d^T^, = , (1) 

for the energy-momentum tensor T^^ = (e + p)u^Ui, — pS^^^, where e and p are the energy 
density and pressure in the comoving frame, and is the 4-velocity of the element of 
the fluid |l2l. The initial state is formed from subcollisions of constituents (nucleons, 
constituent quarks and/or partons) of colliding nuclei and is glue dominated at early stages. 
The LPM effect implies that for a subcollision at the origin, x = {t, z, r) = 0, the secondary 
particle formation vertices lie on a hyperbole of constant proper time r, = — 2;^ ~ Tq , 
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and e,p do not depend on the space-time rapidity rj = ^In (fzf) of the comoving reference 
frame ||^, 0]. In the 1+1-dimensional approximation, this leads to the celebrated Bjorken 
equation 

1^ + ^ = 0. (2) 

OT T 

According to the lattice QCD studies, the familiar |12] = | holds for a velocity of 



sound Cs in the QGP excepting a negligible narrow region of the hadronization transition 
temperature Th ~ IGOMeV and energy density e/j ~ 1.5 GeV/fm^ |14|. With the 
equation of state p = c^e, the Bjorken equation has a solution 



e oc r 



(3) 



Widely varying estimates for emax and the proper time tq are found in the literature [0, |T^, 
|15|, |16[. However, as Bjorken has argued 0, emax oc and the actual dependence of the 
Bjorken lifetime Tb on emax is rather weak: 



tb = To 



n. 



1 



-max \ i+c 



eft 
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For central PbPb collisions, for which there is some experimental evidence for the QGP 
formation 0, the typical estimates are ^ 3f/c at SPS [|14|, |l7l and ^ 6f/c at RHIC 
0, which are much larger than the standard estimate tq ~ 0.5 f/c. 

Now we turn to the major theme of collective transverse expansion, which is driven 
by radial gradient of pressure. As we shall see, the radial flow is nonrelativistic. Then, 
to the first order in radial velocity Ur, the radial projection of (|l]) gives the Euler-Landau 
equation 

/ f)( 1= -i- ( 1= -i- nW fin 

(5) 



dp 
or 



dr ' ' \ dr t ^ 

in which we can use the Bjorken's solution for e and p = c^e. Then the Euler-Landau 

equation can be cast in a simple form 



dur cl cl dlogp 

—Ur = 



(6) 



dr t"'' {1 + c^) dr 
The important point is that the transverse expansion of the QGP fireball can be neglected 
which we can justify a posteriori. For this reason the time dependence of the logarithmic 



derivative D{r,T) = can be neglected, it is completely determined by the initial 

density profile and depends neither on the temperature nor fugacities of quarks and gluons, 
which substantially reduces the model-dependence of the transverse velocity. Then the 
solution of (H) subject to the boundary condition Uri^To) = is 



Ur{r,T) = ^f^dtt~^'^D{r,t) 

J- C„ J TO 



clTD{r,0) 



1 — r4 

TO i C„ 



1- ' 7 



^0 ' ■ 
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The model-independent estimates for the initial density /pressure profile are as yet lack- 
ing. AT RHIC and higher energies of the initial state is expected to be formed by semihard 
parton-parton interactions for which nuclear shadowing effects can be neglected P, p!5| . 



Then for central collisions e(r, Tq) oc T^{r), where k = 2 and T^ir) = J dznj^Sc^ + r^) 
is the density of constituents, Tyi(r) ~ exp(— -^), where Ra ~ l.lA~^/^fm is the nuclear 
radius. In another extreme scenario of strong shadowing and of strong LPM effect the soft 
particle production is not proportional to the multiplicity of collisions of fast partons 0, ^ 
and /c = 1 is more appropriate. Hereafter we take k = 2. In any case, the logarithmic 
pressure gradient is approximately linear, D{r,t) ^ 2kr/R\, and according to the solution 
(0) the displacement of the fluid element Ar is proportional to the radius, Ar oc rr^. 
Consequently, we have the Hubble-type radial rescaling 

\{t)^1 + — ^1+(—)\ (8) 
r \tt/ 

where 

V rvf^s 

has a meaning of the lifetime against transverse expansion. For central PbPb collisions 
eq. (P) gives tt ~ 10k~^'^f/c which is larger than the above cited estimates of tb and at 
SPS and RHIC energies the transverse expansion of the fireball can be neglected. 

In a quasi-uniform plasma the hydrodynamic expansion lasts until t = th- For long-lived 
QGP and r ^ tq we can use the Bjorken's solution 

Si)' 
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which gives the position of the hadronization front 

rhir) = Ra 



l±^ln^ 
k T 



and the radial velocity at the hadronization front 



RA(1-Ci 



4A;(l + c2)ln^ 



r 



i-ci 



(12) 



For the usually discussed tb and tq we have tb ^ tq. For such s long-lived QGP, tb ^ tq, 
the radial velocity takes the maximal value Um at t ~ 1/ ^/e, 



Ur 



'2A;(l + c2 



RAil-ct 
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We 



l-d 
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It is remarkable that the average radial acceleration u^Itb is approximately constant. 
The solid line in fig. 1 shows the maximal velocity Um evaluated from ([T2|). The large-r^ 



approximation ([T3|), shown by the dashed fine, reproduces these results to better than ~ 4% 
at = 3 f/c and better than ~ 1% at = 8 f/c. For the above cited estimates for tb 
in central PbPb collisions we find Um{SPS) ~ 0.13 and Um{RHIC) ~ 0.28, consequently 
the nonrelativistic expansion approximation is justified very well. Now notice, that for a 



long-lived QGP the hadronization front (|TT]), shown in fig. 2, and u{T)/um depend only on 
the scaling variable t = t/tb, with an obvious exception of the short-time region r ~ tq. 
This scaling property is clearly seen in fig. 3 where we show the u/um as a function of 
t = t/tb- Notice a convergence to a universal curve with the increasing tb- 

The most interesting quantity is a radial velocity distribution which can be evaluated 
experimentally from the Doppler modifications of the thermal spectrum. In order to test 
our results one needs particles which are radiated from the hadronization front. The 
standard scenario is that the hadronization transition is followed by an expanding mixed 
phase which, however, does not contribute to the transverse velocity because in the mixed 
phase c^O is negligible small. The mixed phase is followed by a hydrodynamic expansion 
and post-acceleration of strongly interacting pions and baryons until the hadronic freeze- 
out temperature Tj < Th is reached |]18[. This post-acceleration is negligible for weakly 



interacting and 0-mesons, which gives the desired access to the radial fiow at the 



hadronization transition. In the evaluations of the modification of the the thermal spectrum 
and one needs to know the w-distribution weighted with the particle multiplicity. For the 
of constant hadronization temperature contribution of the hadronization surface to 
the particle multiplicity is 

dw oc rh{r)dT . (14) 

Making use of the solutions ( [TT| ) and (|T^, it can readily be transformed into dw/du oc 
rhidu/dr)^-^ . The important point is that because of the above discussed scaling properties 
of the transverse flow, the velocity distribution is a scaling function of x = u/um'- 

dw ^ 1 fix,TB) ^^^^ 

du Um \J\ — X 

where for a long-lived QGP /(x, r^) does not depend on r^. The square-root singularity at 
a; = 1 is a trivial consequence of the vanishing derivative duir) / dr aX r tb / i/e- In Fig. 4 
we show the scaling function f{x,TB) for tb = 6 f/c. We don't show f{x,TB) for other 
values of tb, because the variations from = 6 f/c to 3 f/c to 9 f/c do not exceed several 
per cent and are confined to a narrow region of x < 0.2. The approximation f{x) = 0.5 is 
good for all the practical purposes. 

In conclusion, we would like to argue that the shape of the velocity distribution is to a 
large extent the model independent one. The generic origin of the square-root singularity 
at X = 1 has already been emphasized, the fact that /(O) 7^ is due to a radiation from 
the surface Vh ~ Ra at early stages of expansion. Above we assumed that hydrodynamic 



expansion continues untill the hadronization transition. Following Pomeranchuk |19] one 
can argue that in the non-uniform plasma the hydrodynamic expansion stops when the 
mean free path 

hnt = , ^ X (16) 
n(rc, T)at 

defined in terms of the transport cross section cTf, is larger that the GQP density varia- 
tion length Z}(r, 0)^^. In the partially equilibrated QGP Imt oc . Then the Pomer- 
anchuk condition gives the temperature Tp at which the hydrodynamic expansion stops, 
Tp oc r / E?}s^. The possibility remains open that at early stages Tp > Th, in which case 



dw oc rh{T){Tp/Th)^dT. This enhanced radiation at early stages at slow radial expansion 
but at higher temperatures Tp would mimic radiation at a lower temperature and higher 
radial velocity. This may result in the apparent depletion of /(O); in order to explore this 
possibility one needs better understanding of Imt near the hadronization transition. 

The NA49 fits to the proton, kaon and pion transverse mass rriT distribution in central 
PbPb collisions at SPS assume identical freeze-out temperature for all particle species 
1^. For positive particles NA49 finds Tf = 140 ± 7 MeV and the transverse velocity 
(u) = 0.41 ± 0.11. However, for the one must take the higher freeze-out temperature 
Tf = Th ^ 160 MeV given by the lattice QCD. Because of the anti-correlation between the 
local temperature Tf and {ut), see Fig. 7 in such a fit with larger Tf to the same rriT 
distribution shall yield smaller {u). 
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Figure captions 

Fig. 1 The maximal velocity of radial expansion Um for for central PbPb collisions as a 
function of the expansion time tb- Shown by the dotted line is the large-r^ formula 

Fig. 2 The time dependence of the hadronization front for central PbPb collisions. 

Fig. 3 The converegence to the scaling behaviour of the time dependence of the relative 
radial velocity u/um for central PbPb collisions. 

Fig. 4 The scaling function /(x, tb) of. Eq. (^) is shown for tb = 6. 
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